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Abstract 


sion line, switch gear, and system level microprocessor 
controls. 


Rocketdyne is under a contract from NASA Lewis Re- 
search Center (LeRC) to design, develop, and build 
the electrical power system for the Space Station. 
This system provides for electric power generation, 
conditioning, storage, and distribution. The initial 
configuration uses photovoltaic power generation. So- 
lar dynamic power generation will be added later. 
The power system control is based on a hierarchical 
architecture to support the requirements of automa- 
tion. 

In the preliminary design and technology development 
phase of the program, various modeling techniques 
and software tools were evaluated for the purpose of 
meeting the Space Station power system modeling re- 
quirements. Rocketdyne and LeRC jointly selected 
the EASY5 simulation software, developed by Boeing 
Computer Services, as a system-level modeling tool. 

This paper describes the application of the selected 
analytical modeling approach to represent the entire 
power system. Typical results of model predictions 
are also summarized. The equipment modeled in- 
cludes solar arrays, dc to dc converters, resonant in- 
verters, battery storage system, alternator, transmis- 


During the advanced development phase of this pro- 
gram, several models have been developed using this 
approach. 

Introduction 

Preliminary design of the Space Station Electrical 
Power System (EPS) consists of a hybrid power gen- 
eration system, Solar Dynamic (SD) and Photovoltaic 
(PV), together with an autonomous high frequency 
(20 kHz) Power Management and Distribution 
(PMAD) system. Figure 1 shows the EPS in a simpli- 
fied schematic diagram. To develop and verify such a 
system, two different test beds are operating at LeRC. 
Also, an Integrated test bed (ITB), which incorpo- 
rates both the power system and the hierarchical con- 
trol system, is in the planning stage. 

The modeling approach and simulation software were 
selected to support the preliminary design phases of 
the EPS, including test bed development. Design 
tradeoffs in equipment, controls, or even the entire 
system are supported by analytical models and com- 
puter simulations to minimize development risks. 

This paper describes analytical modeling require- 
ments, summarizes the status of model development, 



Figure 1, Space Station Electrical Power System Simplified Schematic Diagram 





and reviews representative simulation results. Prelimi- 
nary design of the EPS is supported by the develop- 
ment of power system test beds. The test beds will 
also support the assessment of various analytical mod- 
eling techniques, State-of-the-art simulation soft- 
ware, and test validation of selected modeling ap- 
proaches. This paper presents results of this assess- 
ment. Also presented are representative results of the 
modeling effort with special emphasis on control sys- 
tem development requirements and system integration 
analytical support. 

Power System Analytical Modeling Requirements 

Analytical models are needed to support the design 
and development of components and systems. They 
are also essential for the development of controls and 
related software. The type of model, the model com- 
plexity, and the simulation tools utilized depend con- 
siderably on the purposes of the simulation. 

System analysis, modeling of power components, and 
modeling of modules within the system model are 
based on the hybrid multiport representation. Pa- 
rameters for such representation can be a function in- 
cluding nonlinearilies. The EASY5 software is ideally 
suited for this type of representation. 

The terminal characteristics of each function can be 
described by the following equation: 
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The elements within equation (1) can be determined 
by utilizing slate variable approach or frequency do- 
main (transfer function) representation. EASY5 will 
perform the necessary time or frequency domain 
simulations. 

The state variable approach is ideally suited for time 
domain simulations. In it, a linear, time invariant net- 
work is characterized by two equations of the follow- 
ing forms. 

State equation: 

d/dt |X] = I A] [X] + [B] [U] (2) 

Output equation: 

IY) = [CJ|X] + [D] [U] (3) 

Where, 

U = m x 1 vector representing m inputs 
Y = p x 1 vector representing p outputs 


X = n x 1 vector of n independent variables 
A, B, C, D are constant real matrices, called state 
equation matrices, of appropriate dimensions; A is al- 
ways n x n . 

Equation (2) is a set of n first order differential equa- 
tions usually referred to as the state equation in nor- 
mal canonical form. The set of independent vari- 
ables X « (xi, xz xn) is called the state vector. 

The main advantage of the stale variable method is 
the ease with which it can be manipulated by many 
numerical methods of analysis. Furthermore, the 
state variable method can be extended, without much 
difficulty, to deal with nonlinear networks. 

The state equations are solved using numerical inte- 
gration techniques. The choice of a suitable tech- 
nique is influenced by a number of factors, some of 
which can be problem dependent. For example, 
modeling of power electronic devices, when coupled 
with short transmission lines, results in a stiff system 
of circuit equations, i.e., a system with widely sepa- 
rated values of time constants. Also, full order mod- 
eling of power equipment requires handling of highly 
discontinuous functions. A numerical integration 
technique to solve such a system will require a rela- 
tively small value of an upper limit on the integration 
step size in order to ensure numerical stability. 

Two numerical methods of integration were success- 
fully utilized in the power system modeling without 
compromising numerical stability. They are the for- 
ward Euler and the Boeing Modified Stiff Gear. 

Modeling Results 

The EPS component models were developed using the 
hybrid multiport approach described earlier. These 
component models were integrated into subsystem 
level models from the modular elements stored in 
EASY5 Macro Library. The modular approach re- 
duces the time needed for the development of large 
scale system models and permits model test validation 
at subsystem level. Following is a description of repre- 
sentative models and results of simulation. 

Photovoltaic (PV) Source Model 

The baseline PV subsystem topology uses sequential 
shunt regulation during sunlight to maintain array volt- 
age below safe level, provide regulated dc power to 
inverters, and clamp the array at low voltage for safe 
maintenance. Eclipse power is provided from battery 
storage units at regulated level. The combined battery 
charge/discharge unit (BCDU) provides bus regulation 
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during eclipse and individual charge management to 
maximize battery life. 

The Solar cell model[l] is based on the classical equa- 
tion for a solar cell that includes the effect of an 
equivalent shunt resistance. 

I = Ig - !o { exp[q(V + I R*)/( A-k-T)] - 1} 

- [(V + I-Rs) / Rsh] (4 ) 

where: I = current output of the cell 
]g = light generated current 
lo = reverse saturation current 
q = electronic charge 
A = empirical fitting constant (one for 
ideal junction) 
k = Boltzmann constant 
T = absolute temperature, deg.K 
V = voltage appearing at the cell terminals 
Rs = series resistance 
Rsh = shunt resistance 
S = solar insolation 


The equivalent circuit of this model is shown in Figure 



Figure 2. Solar Cell Equivalent Circuit 


Newton iteration method is used to solve for I and V 
for a given load. Effects of illumination and tempera- 
ture on the current and voltage of individual cells can 
be accounted for by correction coefficients. The per- 
formance of the solar array panel can be determined 
from individual cell data using known techniques! 1 ] . 

The solar array power management approach uses a 
pulse width modulated (PWM) DC Switching Unit 
(DCSU) which provides a low dissipation method of 
matching array power delivery lo load demand while 
maintaining voltage regulation (see Figure 3). The 
system consists of the Sequential Shunt Unit (SSU) 
and redundant photovoltaic control elements (PVCE). 
The FVCE provides a bus voltage error signal to the 
SSU’s which shunt excess array power by shorting in- 
dividual solar cell strings to ground. 

In the SSU, each FET switch circuit is driven by a 
separate voltage comparator which compares the com- 
mon bus voltage from a ramp generator. Each ramp 
signal is dc biased (level shifted) to a unique voltage 


level for each comparator. The average current 
through the associated regulating shunt element is pro- 
portional to the PWM duty cycle. The PWM imple- 
mentation simplifies the electromagnetic compatibility 
(EMC) management by synchronization with the dis- 
tribution bus frequency. 



The nonlinear model developed for the SSU is based 
on the state differential equation representation of the 
error amplifier, a functional generation of the 20 kHz 
sawtooth reference signal, and the dc bias to deter- 
mine the appropriate number of shunted solar array- 
strings. 

The computer model implementation of the PV 
Source subsystem is shown in Figure 4. 



Figure 4. EASY5 block diagram of PV Source 
Subsystem Model 


Transient responses of the PV source during startup 
and load switching have been predicted by the model 
including an idealized constant power load represent- 
ing the negative impedance characteristic of a regu- 
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lated inverter load, and by the actual inverter load. 
Details of the inverter model in this simulation are de- 
scribed elsewhere 1 2] 


Simulation resulls shown in Figure 5 represent an opti- 
mum SSU control performance with the inverter input 
current completely filtered. 



Figure 5. PV Source Subsystem 
Predicted Performance 


Modeling of power equipment such as dc to dc con- 
verters, frequency changers, and resonant inverters, 
transmission lines, and remote power controllers 
(RPC) is based upon hybrid multiport representation . 
The parameters for the such representation are devel- 
oped by utilizing the stale-space approach. The ele- 
ments in the hybrid multiport can be nonlinear or time 
varying. 


Generally, two topological methods can be used for 
the formulation of system equations. One such 
method is the constant topology method while the 
other is the variable topology method [3]. In the mod- 
eling of dc-to-dc regulators, frequency changers 
(FCU), and the main inverter unit (MIU), the vari- 
able topology method was used. This method requires 
computation of the state matrices, which have both 
time varying elements as well as dimensions. 


The power stage model of the battery charge convert- 
er (buck regulator) presents three possible topologies 
that exist in course of a complete switching cycle as 
shown in Figure 6. Similar set of equations were de- 
veloped to represent the battery discharge converter 
(boost regulator) power stage. The control circuits 
can be modeled either by transfer functions or by stale 
equations. The state equations were also used to 


model power equipment because of their ability to de- 
scribe the circuit in more detail and also because non- 
linearities such as op-amp saturation and protection 
functions can be easily added. 
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Figure 6. Buck Converter Power State Circuit 
Topology Modes and Model Equations 

Power Management and Distribution (PMAD) 
Model 

The PMAD subsystem will distribute power from the 
sources to the users. The distribution network will 
consist of three major dual power paths. The Power 
Distribution and Control Assembly (PDCA) will serve 
as the primary/secondary distribution tie point. Each 
PDCA will have the ability to isolate load faults, 
switch/balance loads from the dual busses, and can 
shed loads when under an overload condition. The 
widely separated source and load nodes in the system 
will be controlled by a distributed computer system 
with a hierarchical structure. 
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PMAD Subsystem Model 

The PMAD subsystem model incorporates models for 
busses, remote bus isolators (RBI), transmission lines, 
transformers, load converters, switchable user loads, 
and analytical representation of higher level controls. 

The capacitive bus model represents the summation of 
currents, the rate of rise of voltage being determined 
as: dVn / dt = Z i / C. In general, instrumentation 
loads, snubber circuits, and ac filters were added to 
the bus circuit as shown in Figure 7. 



Figure 7. Bus Model Representation 



Figure 8. Generalized Transmission & RBI 

The switching elemenls(SCR’s) within the RBI are an 
integral part of the transmission line model as shown 
in Figure 8. The ac transmission line model (Nominal 
n ) incorporates the necessary elements for the repre- 
sentation of short transmission line transients with se- 
ries R,L and shunt C elements. The switchable im- 
pedance, Rf & Lf, is provided to simulate faults. 

AC Switch Model 

The ac switch consists of solid slate devices and other 
supporting electronics. The ac switch model, however, 
consists of a binary (two valued) series inductance and 
resistance circuit (3J. The two values represent on 
state and off state. The former are chosen to repre- 
sent actual values for the switch itself as well as those 
due to connecting hardware. The latter are chosen to 
represent switch turn off characteristics in a simplified 


manner. Snubber circuits are not included in this 
model. Following a turn off command the current is 
allowed to cross zero. At the end of that time step, 
off state values of the inductance and resistance are 
invoked such that the non zero switch current asymp- 
totically goes to zero with a specified time constant. 
The impedance of the off state switch should be large 
enough to permit the allowable level of leakage cur- 
rent. Figure 9 shows the logic of switch operation as 
incorporated in the model. 



Figure 9. AC Switch Logic 


The transformer model consists of an ideal trans- 
former in series with an R-L circuit on the secondary 
side to represent leakage inductance and losses. Rep- 
resentation of magnetizing characteristics and inter- 
winding capacitance leads to an extremely stiff system 
and therefore difficulty with integration. 

A digital model representation of the higher level con- 
trols has been developed and tested in order to sup- 
port control system development. This model incor- 
porates data transmission from sensors (including de- 
lay time), control algorithms, processor cycle time de- 
lay, and communication among the processors (see 
Figure 10). 

Models of load converters, from 20 kHz to dc as well 
as ac, have been developed under a separate modeling 
effort [2]. 


5 













pwt mt 


I SOLAR DYNAMIC 
MOOULE 


90 

I controller 


INTERFACE 


MOTOR 

CONTROLLER 


photovoltaic 


I rn r 


ALPHA 

XMCT _ 


| INTERFACE! 

lW»T , J 

Tl 

POWER 

(management 

PRQCEIfQRj 


V 

ROCA 

PROCESSOR 

MRSA 

PROCESSOR 


■■ 

■■ 

iKsruzmi 

tmunjzx^mi 



(») 


m 


WTERPACE 


t ARRAY REGULATOR 
RIVER TER 

charge /DISCHARGE CONVERTER 


Figure 10. Electrical Power Control and 
Management System 


Figure 11 shows the model of a Power Distribution 
Control Unit (PDCU) including the EASY5 blocks 
representing the overlaid circuit model. The purpose 



Figure 1 1 . PDCU model with two test 
sources and three test loads 


of the PDCU is to control the state of switches (RBI’s 
and RPC’s), either by command or automatically, 
e.g., to clear a fault. Additionally, the PDCU moni- 
tors currents, voltages, and powers associated with 
these switches. Results of a simulation on the PDCU 
model are shown in Figure 12. The three load cur- 
rents are shown as a function of manual operations of 
the three RPC’s. For reference, one of the source 
voltages is also shown. Two of the RBI’s are also 
manually switched which affects two of the load cur- 
rents as can been seen in the figure. 

Solar Dynamic (SD) Source Model 




Figure 12. PDCU Simulation Results 

poses [4], This alternator was initially selected as the 
SD source for the EPS. In the PMAD System lest 
bed, the SD source was replaced by an SD Simulator 
which consists of BRU-F Alternator powered by an 
air turbine. Since the air turbine lacks any speed gov- 
erning mechanism, speed regulation is achieved by 
means of a balanced, 3 phase, parasitic load control 
scheme. Similar scheme will be used in the actual 
space station. The alternator is equipped with sell ex- 
cited shunt and series field windings which are con- 
trolled by a voltage regulator. The shunt field gener- 
ates nominal terminal voltage while the series field 
compensates for voltage drop due to loading. Figure 
13 is a schematic representation of the SD Simulator. 
The alternator is a 4 pole, 36000 RPM machine, rated 
at 14.3 kVA , 0.75 lagging pf, 120/208 volts. Me- 
chanically, the machine is of the round rotor type but 
magnetically it is of the salient pole type. 


VOLTAGE 

REGULATOR 



Figure 13. Solar Dynamic Simulator schematic 
Alternator Modeling Equations 


LeRC had developed a Brayton cycle Rotating unit 
with foil bearings (BRU-F), turbine driven, 1200 Hz, 
modified Lundell Alternator for space power pur- 


Modeling equations for the alternator as initially de- 
veloped ((5j) are based on the well known d-q axis 
transformation principle which is routinely applied to 
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the analysis of synchronous machines by the utility in- 
dustry. With known information about the circuit 
constants of the stator, rotor, and the damper wind- 
ings of the alternator, equations are developed which 
relate the winding currents, their respective flux link- 
ages, and the voltages at their terminals. 

The EASY5 implementation of the SD source includes 
dynamics of the rotating equipment, voltage regulator, 
parasitic load speed control, and user load [6J. 


will be an essential tool for the support of full scale 
development of the EPS. 

The computing power required for this type of repre- 
sentation is substantial. Therefore, a proper mix of 
detailed models and reduced order models will be re- 
quired for effective and efficient simulation support of 
the program. 
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Conclusions 

The paper presents an overview of the modeling ap- 
proach for support of the Space Station EPS develop- 
ment. Representative simulation results are also pre- 
sented. Test validation efforts, in progress, show that 
the models generally give correct results. The models 
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